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Melting of two dimensional (2D) clusters of classical parti- 
cles is studied using Brownian dynamics and Langevin molec- 
ular dynamics simulations. The particles are confined by a cir- 
cular hard wall or a parabolic external potential and interact 
through a dipole or a screened Coulomb potential. We found 
that with decreasing strength of the inter-particle interaction 
clusters with short-range inter-particle interaction which are 
confined by a hard wall exhibit a re-entrant behavior in its 
orientational order. 



The structural and dynamical properties of small clas- 
sical two-dimensional (2D) clusters have been the sub- 
ject of recent experimental studies and Monte-Carlo 
and molecular dynamics simulations [ i - pXfl . It was found 
earlier that the particles are arranged in shells and that 
melting of finite clusters is a two step process. With in- 
creasing temperature inter-shell motion develops and the 
system losses angular order. Consecutively radial diffu- 
sion switches on and destroys the shell structure of the 
cluster. The spectrum of such clusters was obtained in 
Refs. for different number of particles. The derived 

minimal frequencies and the corresponding energy barri- 
ers showed that 'non-close packed' clusters are unstable 
against inter-shell rotation. High symmetry clusters (i.e. 
the so-called magic number clusters) have energy barriers 
for intershell motion which are several orders of magni- 
tude larger than those for 'non-close packed' clusters. 

Recently, Bubeck et al || observed re-entrant melt- 
ing in two dimensional (2D) colloidal clusters. The clus- 
ters consist of paramagnetic colloidal spheres which were 
confined in a circular hard wall vessel. The external 
magnetic field induces a magnetic moment M in the 
particles and they interact through a dipole potential 
V{fi,fj) = [jqM /4irrh, where /iq is the magnetic per- 
mitivity, and the interparticle distance. The cou- 
pling parameter, which is the inter-particle interaction 
energy measured in units of the particle kinetic energy 
T = V/ksT, characterizes the order of the system. It de- 
creases by lowering the external magnetic field. In Ref. 
H it was found that with decreasing T, first inter-shell 
rotation appears which destroys the angular order of the 
cluster. Further decreasing the parameter T the system 
unexpectantly regained angular order within a narrow 
range of T and then melts when T is further decreased. 
It was suggested that the observed re-entrant melting be- 



havior is due to the increasing role of the radial particle 
fluctuations which is very similar to an earlier investi- 
gation of laser induced melting of 2D colloidal crystals 
IB; 

Earlier theoretical work |5| on parabolic confined clus- 
ters did not find such a re-entrant behavior which sug- 
gests that the shape of the confinement potential may be 
very important. Another difference is that in the exper- 
imental system the particle motion is strongly damped 
because the colloidal particles move in water. In the 
present paper we investigate the mechanism for the re- 
entrant behavior and address the specific role played by 
the type of confinement (i.e. hard wall versus parabolic) 
and of the functional form of the inter-particle interac- 
tion (short range versus long range) on the melting of the 
clusters. 

In our model the particles are confined by a circular 
hard wall potential (V p = for r < R and V p = oo 
at r > R) or by a parabolic potential V p = ar 2 . The 
particles interact through a dipole potential V(fi,fj) — 
q 2 1 | n — fj | 3 , where q 2 — ii M 2 /4tt, or through a 



screened Coulomb potential V(f 



(q 2 / 



) exp(— k\ fi—fj |), with q the 'particle charge', ?i is the 
coordinate of the i th particle, and 1/k is the screening 
length where k = for a Coulomb cluster and we took 
k = 2/clq for the screened Coulomb cluster, where do 
is the mean inter-particle distance. For given type of 
inter-particle interaction and external confinement, only 
two parameters characterize the order of the system: the 
number of particles N and the coupling parameter T. We 
define the characteristic energy of inter-particle interac- 
tion for dipole clusters as E — q 2 /oq and E = q 2 /ao 
for screened Coulomb clusters, where ao = 2R/N 1 / 2 for 
the hard wall and ao = g 2 / 5 a _1 / 5 for parabolic con- 
finement. In the present calculation we define the cou- 
pling parameter as V — q 2 /agfcsT for dipole clusters and 
r = (q 2 /aottBT) cxp(— koq) for screened Coulomb clus- 
ters. In U a different dimensionless parameter T was 
introduced, where V was taken to be the sum over all 
pairs of particles. Our coupling parameter T is a factor 
2.2447 smaller than the one of jgj for N = 29. 

The ratio of the particle velocity relaxation time versus 
the particle position relaxation time is very small due to 
the viscosity of water and therefore the motion of the 
particles is diffusive. In our simulations we will neglect 
hydrodynamic interactions. Following JD| we rewrite the 
stochastic Langevin equations of motion for the position 
of the particles as the equations of motion for Brownian 
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particles 

dfj _ k B T y^ dV(fi,rj) | dV p (fi) | F L 

dt Dam ^— ' c?r dr m ' 

3=1 

where _D is the self-diffusion coefficient, m the parti- 
cle mass, and Fl the randomly fluctuating force act- 
ing on the particles due to the surrounding media. In 
the numerical solution of Eq. (1) we took a time step 
At < 1CT 4 /(™A>), where n = N/(nR 2 ) is the particle 
density and At was varied within a range (0. 02-^0. 05)sec. 
The radius of the circular vessel R = 36fim and the self- 
diffusion coefficient Dq = 0.35/jto 2 / s are taken from the 
experiment || . Following Ref . j| we consider dipole clus- 
ters consisting of N = 29, 30, 34 particles which have dif- 
ferent types of packing. In the ordered state the sys- 
tems of N = 29, 30 particles are arranged in a triangular 
'closed packed' structure having, respectively, the shell 
structure (3:9:17) and (3:9:18) and ground-state energy 
E = 2.2447^o and E = 2. 2798£' . The cluster with 
N = 34 particles (4:11:19) E = 2.4198£ , has a non-close 
packed structure. 

We find first the ground-state configuration using the 
Monte-Carlo (MC) technique. Our results for the mini- 
mal energy configuration coincides exactly with the en- 
ergy found in S for Coulomb clusters and in for 
dipole clusters with parabolic confinement. In the exper- 
iment (^] the system was first equilibrated and after that 
the particle trajectories were recorded during 30 min- 
utes. In the simulation we equilibrated the system for 
about (5 • 10 5 -J- 10 6 ) MC steps after which we started 
with the statistical averaging over time of the different 
observables. To obtain reliable results with small sta- 
tistical error we follow the particle trajectories typically 
during 10 7 time steps. 

We calculated the time dependence of the mean an- 
gular displacement of the particles in a specific shell 

N 

= ~ 0i(O))/W, where 9 t {t) is the angular 

z=l 

position of the i th particle and N is the number of par- 
ticles in the shell. In Fig. 1 the angular displacements of 
the particles in the first shell (the most inner) and the 
second shell are given as function of time for the cluster 
of N — 29 particles. The angular motion in the third 
shell is very small because its motion is hindered by the 
hard wall. Notice that both the first (thick curve) and 
the second (thin curve) shells take part in the angular 
motion, but the former rotation is more prominent. The 
inter-shell motion has no preferential direction and with 
time it can be either a clock wise or a counter clock wise 
rotation. With decreasing coupling parameter T inter- 
shell rotation becomes more pronounced. 

In order to characterize the order of the system we 
calculate the angular diffusion of the particles over a 
30mm x 1000 time interval. The angular diffusion co- 
efficient can be written as 

Dg = (< A6(t) 2 >-< A0(t) > 2 )A, (2) 



where <> refers to a time averaging, and the mean 
relative angular displacement rotation of the first shell 
(0 1 (t)) relative to the second (6 2 (t)) one is defined as 
A6(t) = O 1 ^) - 9 2 (t). The radial diffusion coefficient is 

1 N 

AR2 = a? E(< r ^ « 2 > - < »•<(*) >2 )/ fl o> (3) 

i=l 

which is a measure of the radial order in the system. 

In Fig. 2 the angular and radial diffusion coefficients 
are shown for three different dipole clusters with hard 
wall confinement subjected to the same conditions as in 
the experiment H . For the cluster with N = 29 particles 
the angular diffusion (solid dots in Fig. 2(a)) monoton- 
ically increases with decreasing coupling parameter up 
to r ^ 30 which is a manifestation of angular melting. 
In the interval T = 10 -j- 30 the inter-shell diffusion re- 
mains practically constant, and with further decreasing 
r it is reduced to about a 20% smaller value. In the 
latter region the radial diffusion coefficient starts to rise 
(open dots in Fig. 2(a)), but the cluster retains its shell 
structure. In the range 3 < T < 8 the cluster oscillates 
between the ground state (3:9:17) and the metastable 
state (4:8:17) which leads to a reduction of the angular 
fluctuations. Further decreasing the coupling to T « 5 
both Dg and AR 2 rises quickly, indicating the onset of 
melting. A similar qualitative behavior was observed for 
the dipole cluster with N — 30 particles (see Fig. 2(b)). 
In the non-close-packed cluster with N = 34 intershell 
rotation occurs over all T's considered in the experiment 
(see Fig. 2(c)) and no clear regaining of angular order is 
found in the region 3 < T < 8. 

In order to obtain further insight into this re-entrant 
behavior we investigated the conditions under which this 
novel effect can be observed. Therefore, we varied the 
following parameters: 1) the viscosity of the medium the 
particles are moving in; 2) the way of decreasing the cou- 
pling of the system (i.e. temperature change versus inter- 
particle interaction strength change); 3) the range of the 
inter-particle interaction; and 4) the form of the confine- 
ment potential. 

To study the melting behavior of the cluster under con- 
dition of low viscosity with hard wall confinement we per- 
formed Langevin molecular dynamics (MD) simulations. 
The Langevin equations for the system of -/V interacting 
particles are 

d 2 fj _ _ i/ dn L + F± + Fl_ ^ 
dt 2 dt m m 

where v = Dom/kBT is the friction and E the force 
from the inter-particle interaction. As an example, we 
consider N = 30 dipole particles moving in a medium 
with a viscosity which is 10 4 times smaller than the one of 
water. Such a low viscosity corresponds to the situation 
of colloidal particles moving in a gas with pressure lPa. 
In Fig. 3 the angular diffusion coefficient (solid circles) is 
plotted as function of T. Note that now Dg is about 
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a factor 10 4 larger as compared to previous case (see 
Fig. 2(b)). It is clear that changing viscosity does not 
destroy the re-entrant like behavior. In fact changing 
viscosity will not alter the statistical properties of the 
system but will only change the time scale for relaxation 
to equilibrium. 

The coupling (r) of the system can be varied in two 
different ways. One can decrease the inter-particle inter- 
action as was done in the experiment || or increase the 
temperature of the system in order to induce melting. 
The latter approach was used for the N = 30 cluster and 
the angular diffusion constant is shown in Fig. 3 (open 
circles). The inter-particle interaction was chosen equal 
to the previous one at T = 300 for T — 300K. By heat- 
ing the system no re-entrant melting behavior is observed 
and Dg increases monotonically with increasing temper- 
ature T ~ 1/r. Thus re-entrant behavior is only found 
when the coupling of the system is decreased by chang- 
ing the inter-particle interaction strength at fixed T. The 
reason for this unexpected behavior is that by increasing 
temperature the self-diffusion constant, D$ — ksT/um, 
increases which, because Dg ~ Do, implies that the solid 
dot results in Fig. 3 have to be multiplied with T ~ 1/r, 
and the resulting curve (open dots) does not exhibit a 
local minimum. The numerical results for the dimen- 
sionless variables (e.g. Dg/Dg) depend only on T (and 
N, confinement and inter-particle interaction) but when 
converting them into physical variables T, v, m enters. 

Next we investigated whether the type of inter-particle 
interaction influences the occurrence of re-entrant behav- 
ior. We consider a cluster with long range Coulomb 
interaction (N = 37 (3:9:25) having the same internal 
structure as the N = 30 dipole cluster) with hard wall 
confinement using Brownian dynamics (1). In Fig. 4(a) 
the angular and radial diffusion coefficients are shown as 
function of T. Notice that the Coulomb cluster shows 
completely different melting behavior and Dg increases 
monotonically with decreasing T. Melting takes place at 
r ~ 40. Next we consider short range inter-particle inter- 
action and as an example we took the screened Coulomb 
cluster (k — 2/ao) with N = 30 (3:9:18) particles. In 
Fig. 4(b) the angular and radial diffusion is shown as 
function of the coupling. Dg exhibits a clear re-entrant 
behavior which correlates with an increase of the radial 
diffusion. Thus we may conclude that only clusters with 
short range inter-particle interaction in a hard wall vessel 
shows re-entrant behavior. 

Last we study the effect of the shape of the confinement 
potential and consider the melting behavior of a cluster 
with short range inter-particle interaction in a parabolic 
well. We choose the N = 25 dipole cluster (3:9:13). Fig. 5 
shows the angular diffusion and the radius of the cluster 
as function of T. Dg clearly does not exhibit a local min- 
imum, it rises uniformly with decreasing T, and melting 
occurs for T « 5, as in the case of hard wall confinement. 
The radius of the cluster R and the mean inter-particle 
distance changes proportionally to T 1 / 2 . 

In conclusion, we studied the melting transition of 2D 



clusters with dipole and screened Coulomb type of inter- 
action confined by a hard wall or a parabolic external po- 
tential using Brownian dynamics simulations. Langevin 
molecular dynamics simulations were carried out in the 
case of small viscosity. We found that only clusters with 
short range inter-particle interaction and confined by a 
hard wall well exhibit angular freezing before melting, ir- 
respective of the value of viscosity. But it is essential that 
the coupling of the system is reduced by decreasing the 
inter-particle interaction at constant temperature as was 
done in the experiment ||. In the other cases, either of 
Coulomb clusters or with parabolic confinement the sys- 
tem shows the usual two step melting behavior without 
any re-entrance. 

We showed that re-entrant behavior is a consequence 
of the interplay between angular order and radial os- 
cillations where an increase of the radial fluctuations is 
able to induce angular order in clusters with magic num- 
ber. With decreasing T, first angular motion sets in, 
because it is governed by the lowest energy barriers. Fur- 
ther increasing T leads to an increase of the radial mo- 
tion/fluctuations which hinder the angular motion. The 
later prevents angular motion in case of hard wall confine- 
ment. But for parabolic confinement the average inter- 
particle distance (see Fig. 5) decreases which results in 
an equal scaling of the energy barriers for inter-shell and 
intra-shell motion. This contrasts with the hard wall 
confinement case where the inter-particle distances are 
unaltered and the energy barrier for inter-shell jumps 
decreases leading to an increase of the radial fluctua- 
tions and the inter-shell jumps. Thus anharmonic effects 
are essential for the occurrence of this re-entrant behav- 
ior. Anharmonicity is enhanced in systems with hard 
wall confinement and for short-range inter-particle inter- 
action. 
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FIGURES 



FIG. 1. The angular displacement of the particles on 
the first (thick curves) and second (thin curves) shell for 
a cluster with N — 29 particles for different values of the 
coupling parameter: (a) T = 150, (b) T = 46, (c) T = 10, 
and (d) r = 3. 

FIG. 2. The angular diffusion Dg (solid circles) and 
the radial diffusion (Ai? 2 ) (open circles) coefficients as 
function of T for clusters with (a) N = 29, (b) N = 30, 
and (c) N — 34 particles. 

FIG. 3. The angular diffusion coefficients as a func- 
tion of the strength of the inter-particle interaction at 
constant temperature (solid circles) and as function of 
temperature (r ~ 1 /T) for fixed inter-particle interaction 
strength (open circles) for a dipole cluster with N = 30 
in the case of small viscosity. 

FIG. 4. The angular (solid circles) and radial diffu- 
sion (open circles) coefficients as function of T for: (a) 
the Coulomb cluster with N — 37 particles, and (b) the 
screened Coulomb cluster with N = 30 particles confined 
by a hard wall potential. 

FIG. 5. The angular diffusion coefficient (solid circles) 
and radius of the cluster (open triangles) as function of 
r for a N = 25 dipole cluster in a parabolic well. 
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